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Abstract Radio observations of solar flares often reveal various periodic or 
quasi-periodic oscillations. Most likely, these oscillations are caused by magneto¬ 
hydrodynamic (MHD) oscillations of flaring loops which modulate the emission. 
Interpretation of the observations requires comparing them with simulations. 
We simulate the gyrosynchrotron radio emission from a semi-circular (toroidal¬ 
shaped) magnetic loop containing sausage-mode MHD oscillations. The aim is 
to detect the observable signatures specific to the considered MHD mode and 
to study their dependence on the various source parameters. The MHD waves 
are simulated using a linear three-dimensional model of a magnetized plasma 
cylinder; both standing and propagating waves are considered. The curved loop is 
formed by replicating the MHD solutions along the plasma cylinder and bending 
the cylinder; this model allows us to study the effect of varying the viewing angle 
along the loop. The radio emission is simulated using a three-dimensional model 
and its spatial and temporal variations are analyzed. We consider several loop 
orientations and different parameters of the magnetic field, plasma, and ener¬ 
getic electrons in the loop. In the model with low plasma density, the intensity 
oscillations at all frequencies are synchronous (with the exception of a narrow 
spectral region below the spectral peak). In the model with high plasma density, 
the emission at low frequencies (where the Razin effect is important) oscillates in 
anti-phase with the emissions at higher frequencies. The oscillations at high and 
low frequencies are more pronounced in different parts of the loop (depending 
on the loop orientation). The layers where the line-of-sight component of the 
magnetic field changes sign can produce additional peculiarities in the oscillation 
patterns. 
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1. Introduction 


Waves and oscillations are a common phenomeno in solar flares. They are be¬ 
lieved to reveal themselves as quasi-periodic pulsations in spatially unresolved 
observations.^tdifferen^^avelengths_[see^_ej^jj.^he.^evie^.o^Nakariakov and 
Melnikov. 120091 1 and have b een directly detected us i ng the EUV observations 


with high spatial resolution ( 

Nakariakov et all 19991: Aschwanden et al. 19991: 

De Moortel and Nakariakovl. 

2 OI 2 II. In particular, the observed periodicities of 


the quasi-periodic pulsations in the radio emission, hard X-rays, and 7 -rays (from 
a fraction of a second to several minutes) agree with the expected magnetohydro¬ 
dynamic (MHD) timescales which strongly suggests a close relationship between 
the pulsations and the MHD waves and oscillations ( Nakariakov and Melniko^ . 


2 OO 9 I : De Moortel and Nakariak^ . 1201^ 1. It should be noted, however, that the 


solar quasi-periodic pulsations are usually superimposed on a rapidly varying 
background, so that recovering their parameters is not easy; some previously 
repo rted conclusions may be unreliable due to the methodology used I Gruber 
et ai, 2 OIII I. 

Observations of the waves and oscillations in the flares can be used to diagnose 
the quasi-periodic magnetic reconnection pro cesses that are believed to be th e 
most probable source of these oscillations { e.g., Kliem. Karlickv. and Ben3 . l2000ll . 
The MHD waves and oscillations can also be used to diagnose the coronal plasma 
and magnetic,feld^£arameters_[^koronaI,seismolog^]]j_see 4 _ej£;^Nakariakov and 
Verwichte, 200!tI1 . Moreover, the waves and oscillations are beli~ed to pla y a key 
role in the energy transfer both in the quiet corona {e.g., Erdelvil . I2OO6I I and in 
the flares ([Fletcher and Hudsonl . I2OO8I): they can modula te and/or trigger the 
magneticreconnectioiygrocessesj Chm^n^^ries^j200^Doyle, Popescu, and 
Tarovan. l200mNakariakov et all 200fil l and thus affect greatly the behaviour of 
unstable active regions and development of the solar flares. 

One of the potential observable manifestations of MHD waves is quasi-periodic 


pulsations of radio emiss i on (e. q ., lAschwandenl . 119871 : iNakariakov and Melnikov 


2OO9I : Kupriyanova et ai . 201(TI 1. because these waves can modulate both the co¬ 
herent and incoherent radio emission mechanisms. Interpretation of the observed 
quasi-periodic pulsations requires comparing the observations with simulations 
in order to determine the origin of the pulsations (because they can be caused 
not only by the MHD waves but also, e.g., by quasi-periodic injection of energetic 
particles); if the oscillations are due to waves, it allows us to identify the wave 
mode and other characteristics of the MHD waves. In turn, identification of the 
wave mode is of high importance for the development of the coronal seismology. 
As a rule, the basic characteristics of the pulsations (such as the period) have 
no unique interpretation and therefore some more subtle features (such as the 
relations bet ween the pulsations observed in different spectral channels) have to 
be an alyzed ( Fleishman. Bastian. and Garvi . 20081 : Mossessian and Fleishman . 
2ni2ll . 

Perhaps the best tool to study the physical parameters and processes in solar 
flares is their microwave emission, because it is produced mainly due to the 
incoherent gyrosy nchrotron radiation mechanis m which is (a) relatively simple 
and well studied ( Melrosel . Il968l : Ramatv . IQBOh and (b) highly sensitive to the 
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Microwave emission from an oscillating 3D magnetic loop 


magnetic field and plasma parameters. There are a lot of papers dedicated 
to simulating the gyrosynchro tron emission in solar flares. In par t icular , two- 
dimensiona l 1201 models Ic.g.. fAlissandrakis and P rel^-Papademal.ll984l: Kleiii 
and Trottet.ll984IPreka-PaDadema and AlissandrakiJ 19881 : Kucera et all 1993 : 


Kundu et all 120011) and three-dimensional ISDI models (e.a.. Preka-Pamdema 


and Alissandrakis , llOOd: iBastian. B enz, and Garvl . Il998l : iNindos et oil I 2 OOOI : 


Simoes and Costal . 2000l : Costa et a/.l . l201,'lh have been employed to account for 


the structure of the emitting region. Some simulations involve the magnetic field 
struc ture obtained by extrapolating the observed photospheric mag netograms 


fe.o.. IGarv. Fleishman, and Nital . l2013l : lKuznetsov and Ko ntan. 1201511 . take into 


account anisotropy of the emitting electrons (e. a., FleishmanandMelniko^ 2003 
Altvntsev et al . 2008; Simoes and Costa . 2010l : Kuznetsov. Nita. and FleishmarJ . 
201 il l and/or evolution of the ekctrondistrihutiona_dmmg_^mes_[e;^^Tzatza- 
kis, Nin dos. arid Alissandrakis. 20081: Reznikova et al . 20091: Kuznetsov and 
Zharkova. I 2 OIOI : iKuznetsov and Melnikovl .' 2012 ll . However, possible modulation 
of the emission by MHD waves has been studied only in a few work s and by 
using relatively si mple models (iKoov l ova. Stepanov, and T saol. 20021: Nakarl 


akov and Melni kov, 20061: Reznikova et al. 2007: Mos sessian and Fleishman! . 
2 OI 2 I 1 . Recently, iReznikova. Antolin. and Van Doorsselaereil2014 ) have simulated 
the gyrosynchrotron radio emission from an oscillating magnetic cylinder (with 
standing sausage waves). It was found that even for the simple model considered 
(a straight cylinder), calculating the emission requires 3D simulations that also 
take into account the spatial structure of the MHD waves, because the results can 
differ fromsim 2 le.^stimates_Jor^^_S 2 atiall^Miomo 2 eneous_source_[c^Reznikova 
et al., 2 OO 7 I : Mossessian and Fleishman! . 20121) . At the same time, another com¬ 
plication arises from the fact that the solar flaring lo o ps are (in contrast to the 
model of Reznikova. Antolin. and Van Doorsselaer^ ( 2014 11 evidently curved. 
This work is dedicated to simulating the radio emission from such curved loops. 

Simulation of MHD waves in flaring loops is a nontrivial task. Oscillation 
modes of a s traight overdense cvlindej^ of infinite length are well studied l Edwin 
and Roberts, 19831 : iNakariakov and Verwicht J 20051 . etc.), and the numerical 
simulations can be greatly simplified by using t he analytical results. However , 


these results are not applicable to curved loops ( Van Doorsselaere et a,L . 2004l . 


Variation of the plasma and magnetic field para meters as well as of the loop 
radius withHieight_ 2 resents_^ven^more_diffipultie^ Aiidnes_^t_a|. ^ 200^Andries 


and Gaily, 1201111 . Therefore in this work we use a simplified model in which an 
oscillating cylinder is artificially bent to form a semi-circular loop without actu¬ 
ally using the wave relations in a semi-torus; the magnetic field, plasma density 
and temperature, and the loop radius (in the equilibrium state, i.e., without 
oscillations) are assumed to be constant along the loop. The main aim of this 
model is to study the effect of varying the viewing angle (i.e., the angle between 
the local magnetic field and the line of sight) along the loop. Although such a 
model cannot correctly describe the entire loop, it can reproduce the effects of 
the MHD oscillations in local parts of the loop, where the effects of the loop 
curvature and inhomogeneity on the MHD wave parameters are negligible. We 
consider the incoherent gyrosynchrotron radiation of energetic electrons which 
makes the main contribution to the flare emission in the microwave range. 
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The simulation model is described in Section [2J The simulation results are 
presented and discussed in Section [31 The conclusions are drawn in Section jT] 


2. Model 


2.1. MHD Waves 


We use the model described earlier bv Zaitsevand.Ste2ano^j[W75|^Edwin and 
Roberts (ilQS.li l. Antolin and Van Doorsselaerel (l2013l b etc. In this model, there 
is an overdense cylinder with the magnetic field directed along its axis. The 
plasma density and temperature inside the cylinder are higher than outside, 
while the magnetic field strength inside the cylinder is lower than outside, to 
maintain the pressure balance. We consider the azimuthally symmetric case when 
all parameters depend only on the distance along the cylinder axis z and the 
distance from this axis r. In an equilibrium case, the magnetic field, plasma 
density, and temperature are assumed to be constant inside and outside the 
cylinder; they experience a sharp jump at the cylinder boundary (at r = R). 

The formulae describing the MHD disturbances ar e given , e.g., in the works 
of IZaitsev and Stepanov! (Il975[). Edwin and Robert^ (Il983ll. Antol i n and Van 
Doorsselaere ( 20I3f l. and Reznik ova. Antolin, and Van Do orsselaere ll20I4 b see 
especially the paper of Reznikova, Antolin. and Van Doorsselaere ( 20I4l l where 
the solutions of those equations are presented for the same conditions as used 
in this work. Therefore we describe the properties of the MHD waves only 
briefly. We consider the sausage mode - the compressible azimuthally symmetric 
(to = 0 ) oscillations that affect both the magnetic field and plasma density and 
temperature. We consider both the standing and propagating waves which are 
described by the relations 


Astand oc sin(a;t) sin(/cz + Aprop oc sin(a;t - /cz + (1) 

respectively, where A is the wave amplitude {e.g., in velocity), w and k are the 
wave frequency and wavenumber, respectively, and ip a is the initial phase. 

Eigure[T]presents an example of how the parameters vary with the coordinates 
in a standing wave (the oscillation phase is chosen to provide the maximum 
deviation from equilibrium). In general, the axial magnetic field B^, plasma 
density np, and temperature Tq vary in phase and their variations are maximal 
at the cylinder axis (we are interested only in the parameters inside the cylinder, 
because, as shown below, this is the region where the radio emission is assumed 
to be generated). Variations of the above parameters decrease with increasing r 
and become zero at some radial distance r* nearby the cylinder boundary (r* ~ 
0.9i? for the parameters considered here). Then the amplitude of the variations 
increases slightly, but they are now in anti-phase with the variations of the same 
parameters in the inner part of the cylinder. Variations of the radial magnetic 
field Br are shifted by 1/4 period with respect to other parameters; they are 
zero at the cylinder axis and reach their maximum at r = r* (without the phase 
reversal at that radial distance). Since we consider low-amplitude oscillations. 
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Microwave emission from an oscillating 3D magnetic loop 




Figure 1. Profiles of the magnetic field strength (a) and direction relative to the tube axis 
(b) in the oscillating magnetic tube at different distances from the tube axis. The standing 
wave, low-density model (see Tabled is considered; one wavelength is shown. Profiles of the 
plasma density and temperature are similar to those of the magnetic field. 


we have \Br\ <C \Bz\ (we remind that Bz includes here both the equilibrium 
magnetic field Bq and the perturbed field component i.e., B^ = Bq + SB^)- 
The radial component Br = SBr makes almost no contribution to the resulting 
magnetic field strength; however, it affects slightly the magnetic field direction 
causing it to deviate from the cylinder axis. Variations of the cylinder radius due 
to the MHD oscillations are very small and do not affect the radio emission. 

The propagating sausage wave looks like the structure shown in Figure [1] 
moving along the z axis. All other properties of this wave (including the phase 
relations between different parameters) are the same as for the standing wave. 

We c onsider two sets of model parameters, the same as have been used in the 
work of iReznikova. Antolin. and Van Doorsselaerel ( 2014l l: the parameters are 
listed in Table [TJ The parameters were chosen to provide the same wavelength 
A for the given cylinder radius R. In the calculations below, we use mainly the 
low-density model because those conditions are more favorable for generation 
of the gyrosynchrotron radio emission. In the high-density model, the emission 
intensity is usually lower; at low frequencies, the emission becomes strongly 
affected by the Razin effect (see Section III. 4. 31) . 


2.2. Source Geometry 


As stated before, the aim of this work is to study the effect of the loop curvature 
on the radio emission. For this, we have created model semi-circular (semi- 
toroidal) loops by bending the above-mentioned overdense cylinder (see Figure 
[21). In this case, the longitudinal coordinate within the cylinder z is transformed 
into the distance along the loop axis, and the radial coordinate r - into the 
distance from the loop axis (the explicit formulae are given in the Appendix); 
this allows us to obtain the plasma and magnetic field parameters at any given 
point. Since the lengths of the solar flaring loops are usually much larger than 
their radii, we have created the models by stacking several MHD simulation 
blocks (each corresponding to one wavelength, like in Figured]) along the loop. 
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Table 1. Parameters of the MHD models: a) cylinder radius (R), magnetic fields (B), plasma 
densities (no), and temperatures (To) inside (“in”) and outside (“out”) the cylinder in the 
equilibrium state, concentration of the energetic electrons (nt) inside the cylinder in the 
equilibrium state; b) wavelength (A) and period (P) of the sausage mode, maximum deviations 
of the magnetic field strength (SB), magnetic field direction (SO), plasma density (5no), plasma 
temperature (<5To), and concentration of the energetic electrons (^n^) inside the cylinder from 
the respective equilibrium values. 


Model 

Low-density 

High-density 

R [km] 

Equilibrium cylinder parameters 

1000 

1000 

Bin [G] 

160 

50 

Bout [G] 

161 

56 

no in [cm-3] 

4 X 10® 

lOio 

no out [cm“3] 

10® 

3 X 10® 

To in [K] 

10^ 

10^ 

To out [K] 

2 X 10® 

2 X 10® 

»^bin [cm-3] 

10® 

2.5 X 10® 

A [km] 

Wave parameters 

2825 

2825 

P[s] 

0.35 

1.5 

(5Rin)max [G] 

4.05 

2.78 

((50in )max 

0.69° 

1.52° 

(5?l0in)max [cm ^] 

1.0 X 10® 

6.0 X 10® 

((^Tq in)max [K] 

1.7 X 10® 

4.1 X 10® 

(<inbin)max [cm-^] 

2.5 X 10® 

1.5 X 10'^ 


As the “basic” model, we have used the loop containing three wavelengths; such 
loop has the length of 8475 km and the height {i.e., the curvature radius) of 
about 2700 km. 

The loop orientation is described in general by three Euler angles. However, 
since rotation around the line of sight does not provide us with any new infor¬ 
mation (it results simply in the same rotation of the simulated radio images), it 
is suf5cientJg^jgasideLlh£:^liMig^:gl^alzJ^Q:^Sg=i^^^'^®^®°’^’ Nita, and 
Fleishman, 1201111 . In this work we assume that the loop is vertical and located 
at the solar equator. The loop orientation is described by two angles: the angle ip 
between the loop plane and the equatorial plane and the longitude A. A possible 
loop inclination is not considered for the same reason as described above: for 
a vertical loop, the variation of the angles il> and A is sufficient to obtain any 
possible loop orientation with respect to the observer (with an accuracy up to 
rotation around the line of sight); therefore our results are applicable to inclined 
loops as well. However, we note that the loop locations at the solar disk which 
are discussed below should be treated as approximate (illustrative), because they 
correspond to the zero inclination only. 

It is not possible to present here the results for all possible orientations; 
therefore we have focused on three illustrative cases (see Figure [5] and Tabled]): 
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Disk center Limb A Limb B 



Figure 2. Three orientations of the model magnetic loop considered in the paper: 3D views 
of the loop (top row); profiles of the viewing angle (relative to the magnetic field, in the 
equilibrium state) along the loop axis (middle row); simulated profiles of the radio intensity 
in the equilibrium state at different frequencies (for the low-density model and the electron 
spectral index of <5 = 4) along the loop axis (bottom row). The visible coordinate along the 
loop s is measured from the upper (as shown in the 3D views) footpoint. 


a) The loop is located near the disk center. In this case, the viewing angle 9 
(relative to the magnetic field) varies in a wide range along the loop, being close 
to zero (or 180°) in the footpoints and being equal to 90° at some point near 
the loop top (see Figure [5]). The coordinate s in Figure [5] is the visible distance 
along the loop axis, i.e., it takes into account the projection effects. 

b) The loop is located near the solar limb and its plane is almost parallel 
to the equatorial plane (“Limb A” orientation). In this case, in contrast to the 
previous one, the viewing angle is close to 90° in the loop footpoints and becomes 
very small near the loop top. 

c) The loop is located near the solar limb and its plane is almost perpendicular 
to the line of sight (“Limb B” orientation). In this case the viewing angle is always 
close to 90° and weakly varies along the loop (although it is not constant, which 
makes this model diff erent from those studied in the work of Reznikova, Antolin, 
and Van Doorsselaere ( 20141 11. 
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Table 2. Orientations of the magnetic loop used in the simulations: the angle relative to the 


equatorial plane ('j/’) and the longitude (A). 

Orientation 

Disk center 

Limb A 

Limb B 

i/i 

60° 

20° 

60° 

A 

20° 

80° 

80° 


2.3. Electron Parameters 


We assume that the gyrosynchrotron radio emission is produced by the energetic 
electrons with the power-law distribution over energy, f{E) oc E~^] in the 
calculations, the energy range is taken to be 100 keV < E < 10 MeV which 
is typical for the solar flares. The electrons with lower energies are also present 
in flares, but their contribution to the gyrosynchrotron emis sion is m i nor (m ostly 
in the form of self-absorption at low frequencies; see, e.g., Holman ( 2003h l and 
usually can be ignored. The pitch-angle distribution of the electrons is assumed 
to be is otropic (ie., we assume that it is not affected sign ificantly by the MHD 
waves). iReznikova. Antolin. and Van Doorsselaer^ ( 2014l l have considered the 
case when the concentration of the energetic electrons rib inside the cylinder is 
proportional to the thermal plasma density, ribin/f^oin = const. Alternatively, 
one may consider the model in which the energetic electrons are “frozen” into the 
magnetic field, so that their concentration is proportional to the field strength, 
nh in/Sin = const. However, since in the sausage waves the dominant longitudinal 
component of the magnetic field and the plasma density ng oscillate in phase, 
both approaches provide almost identical results. Therefore in the calculations 
to follow we assume that the concentration of the energetic electrons in the loop 
is always proportional to the thermal plasma density and nbin/?^oin = 0.025; 
in particular, in the equilibrium state we have rib in = 10® and 2.5 x 10® cm“® 
for the low- and high-density models, respectively. We assume that there are no 
energetic electrons outside the loop (rib out = 0). 

The gyrosync hrotron radio emission is calcula ted using the Fast Gyrosyn- 
chro tron Codes dP leishman and Kuznetso\ ], |2010|; Kuznetsov. Nita, and Fleish¬ 
man, 


2011: Nita et al 


201 ^; 


these codes include the contribution of the thermal 
free-free emission as well. We neglect the possible harmonic structure of the 
gyrosynchrotron emission at low frequencies by using the “continuous” regime 
of the fast gyrosynchrotron codes, because in real observations this structure 
should be usualH_smoothed__b^_^he_source_inhomo2eneit^_(Kuznetsov, Nita, 
and Fleishman, 120111 1: using the “continuous” approximation also speeds up 
the calculations appreciably. The codes calculate the emission characteristics 
by numerical integration of the radiative transfer equations for the ordinary 
and extraordinary electromagnetic modes while taking into account the possi¬ 


ble mode-coupling effects, 
transverse magnetic field 


i.e., in t eraction between the modes in the reg ions of 
( Cohenl . llObOl : IZheleznvakov and Zlotnifl . 196^ . This 


approach is sufficient to reproduce accurately the emission intensity in both 
the quasi-longitudinal and quasi-transverse (with respect to the magnetic field) 
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Frequency [GHz] 


Figure 3. Total (spatially integrated) radio spectra in the equilibrium state. Solid lines: 
low-density model; dashed lines: high-density model. Two electron spectral indices (5 = 4 
and 6) and the “Limb B” orientation are considered. 


propaKa tion regime s as we ll as the circular polarization in the quasi-longitudinal 
regime ( Nita et ai . 2015h . although simulating the polarization in the regions 
with 9 ~ 90° may req uire considering the full set of Stokes parame t ers (s ee 
the upcoming paper of Reznikova, Van Doorsselaere. and Kuznetsov! ( 2015l ll. 
Only the emission and absorption processes inside the curved magnetic loop are 
considered, i.e., the emission outside the loop is assumed to propagate like in 
vacuum. 


3. Results 

3.1. Equilibrium State 

Firstly, we present a brief overview of the emission produced by the loop in 
the equilibrium state (be., without MHD oscillations). This emission is very 
similar to the emission from an oscillating loop averaged over the wave period. 
Figure [3] demonstrates the total (spatially integrated) intensity spectra for the 
loop located at the limb (“Limb B” orientation). For other orientations, the 
intensity is slightly lower, but the spectra shapes are qualitatively the same. 
The concentration of the energetic electrons was chosen to provide that the 
spectral maximum (for the low-density model and the spectral index of the 
energetic electrons of d = 4) is located at about 10 GHz - a typical value for the 
continuum microwave emission of solar flares. 

The bottom row in Figure [3] demonstrates the simulated distributions of the 
gyrosynchrotron radio intensity in the equilibrium state along the loop (for the 
low-density model). The intensity distributions depend on the emission frequency 
and loop orientation: for the loop located near the disk center, the optically thin 
emission {e.g., at ft! 25 GHz) is sharply peaked near the loop top; the same effect 
is seen for the emission near the spectral peak (at ft 10 GHz), while the optically 
thick emission {e.g., at ft 4 GHz) is weakly peaked at the loop footpoints. For 
the “Limb A” orientation, the situation is opposite: now the optically thick 
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emission is peaked near the loop top, while the optically thick emission and the 
emission near the spectral peak are peaked at the loop footpoints. A similar 
pattern, although with weaker variations, can be seen in the “Limb B” case. 
In all cases, there are small narrow peaks/dips at the points where the local 
viewing angle equals 90°; they are caused by different polarization characteristics 
of the ordinary and extraordinary electromagnetic modes and are typical of the 
gyrosynchrotron radiation. 

3.2. Intensity Modulation: 2D Images 

Oscillation-induced variations of the radio emission in different parts of the 
magnetic loop tend to compensate each other, thus reducing the variations of the 
total (spatially unresolved) emission. Due to the variation of the viewing angle 
(and hence of the local emission parameters) along the loop, the mentioned 
compensation is not complete. Nevertheless, the variations of the total intensity 
are always much weaker than those of the local (spatially resolved) emission; 
amplitude of the total intensity variations decreases with a decrease of the 
relative wavelength of the MHD waves (le., with an increase of the number 
of the waves in the loop). Therefore we focus below on the characteristics of the 
spatially-resolved emission. 

Figure m demonstrates the simulated radio images at the frequency of 25 GHz 
(optically thin emission) for the standing wave. Two opposite oscillation phases 
are shown. To make the effect of the MHD waves more pronounced, the intensity 
variations were enhanced as 


/^(/) + 5 ( 7 -(/)), 

where (7) is the intensity averaged over the wave period P: 


( 2 ) 



(3) 


We can see in the figure that for the “Limb B” orientation (with an almost 
constant viewing angle) the visible emission structure, in general, reflects the 
structure of the MHD wave: there are three intensity maxima corresponding 
to the points of maximum magnetic field (and hence the maximum concen¬ 
tration of energetic electrons). For other orientations, due to the projection 
effects, the intensity maxima/minima slightly deviate from the magnetic field 
maxima/minima; also, not all MHD waves are visible (*.e., both the intensity and 
the oscillation-induced variation amplitude depend strongly on the coordinate 
along the loop). For the propagating MHD waves, the radio images are similar 
to those shown in Figure SI but the brighter/darker spots drift gradually along 
the loop with time. 

3.3. Phase Relations and Modulation Depth 

We are interested primarily in the qualitative effects that could reveal the basic 
MHD wave characteristics (such as the wave mode). As mentioned above, in 
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Figure 4. Two-dimensional images of the microwave emission (at the frequency of 25 GHz) 
for the standing wave. The low-density model is used, the loop contains three sausage waves 
and the electron spectral index is <5 = 4. Two oscillation phases corresponding to the maximum 
deviations from the equilibrium are shown. The asterisks (*) and diamonds (0) mark the points 
with the maximum / minimum magnetic field at the loop axis. 


real flaring loops the MHD wave properties (including the phase speed) can 
vary strongly along the loop, which is not considered in our model. Therefore it 
makes no sense to compare, e.5., the simulated time profiles of the emission at the 
loop top and the footpoints. Instead we consider the local emission parameters, 
i.e., the temporal delays (phase differences) between the emissions at different 
frequencies at a given point. 

Figures [5]l6] show the time-distance plots of the normalized intensity variation 
(I —(/))/ (/) at the axis of the magnetic loop {i.e., along the line formed by the 
projection of the loop axis on the image plane) for the low-density model. Four 
frequencies are shown, representing the optically thick case (4 GHz), “moder¬ 
ately” optically thick case (6.3 GHz), spectral peak (10 GHz), and optically thin 
case (25 GHz). Above each panel, the respective maximum modulation depth A 
is given; this value is defined as 


A = max 


I-(I) 
(!) 


(4) 
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Figure 5. Time-distance plots of the intensity variation: the normalized intensity variation 
(/ — (/))/ (!) at the loop axis vs. the visible coordinate along the loop axis and the oscillation 
phase. The low-density model (standing wave) is used; the loop contains three sausage waves 
and the electron spectral index is 5 = 4. For each emission frequency and loop orientation, the 
maximum modulation depth A is given above the corresponding panel. 


Note that the maximum modulation depth A is different from the average mod- 
ulation depth e given by Equation ([5]) (which will be used in Section I3.4.2|) . 
According to Figures [S]ini the maximum modulation depth A varies from < 1% 
at low frequencies to > 5% at high frequencies. These values are relatively low 
and, in fact, it would not be easy to detect such oscillations in real observations. 
However, this is because we consider the relatively low-amplitude MHD waves: 
the variations of the emission source parameters SB/ (B) and Srih/ (n-b) do not 
exceed 2.5% and thus the variations of the intensity are comparable to the 
MHD wave amplitude. The qnasi-periodic pulsations reporte d by Nakariakov. 
Melni kov, and Reznikova ( 2003l l. Melnikov et all (|2005l l. and iReznikova et all 
( 2007l l (with the modulation depths of up to « 0.3) require much more powerful 
MHD waves; nevertheless, if those pulsations are caused by sausage MHD waves, 
their characteristics should be qualitatively similar to the case studied in this 
work. Note also that the intensity variations are more pronounced in the optically 
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Figure 6. Same as in Figure [5] for the propagating wave. 


thin fr eque ncy range, in agreement with the results ofiMossessian and Fleishman 
( 2012f ] and Reznikova. Antolin. and Van Doorsselaer j ( 2014l] . 

Figure [S] demonstrates the simulation results for the standing wave. We can 
see that the relative modulation depth, in general (excluding the regions with 
9 ~ 90°), follows the same trends as the unperturbed intensity (see Figure 
HD: for the loop near the disk center, the modulation depth increases towards 
the loop top for the optically thin emission (> 10 GHz) and increases towards 
the loop footpoints for the optically thick emission (< 10 GHz); for the loops 
near the limb, the trends are opposite. Oscillations at low (< 4 GHz, the first 
row of the figure) and high (> 10 GHz, third and fourth rows) frequencies are 
synchronous. The fact that t he optically thin and thick emissions osci l late i n 
phase agrees with the results of Reznikova. Ant olin . and Van Doorsselaerel ( 2014 ), 
but contradicts the c onclusions of | Mo ssessian and Fle ishman ( 2012l ]. The reason 
is that, in contrast to Mossessian and FleishmanI ()2012[ ). we take into account the 
modulation of the viewing angle by the sausage wave but neglect the modulation 
of the energy of the nonthermal electrons. In addition, we use the spatially 
inhomogeneous model in which the source parameters vary along the line of sight; 
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this is especially important for the optically thick emission. S ince two different 
approaches (ours and that of Mossessian and Fleishmiml ( 2012l l'l provide opposite 
results, we conclude that the sausage-mode MHD waves (at least, in application 
to the radio intensity modulation) cannot be treated without a spatially resolved 
3D model. Note that this behaviour (in-phase variations of the optically thin and 
thick emissions) is observed for all loop orientations and at all locations along 
the loop; i.e., it is independent of the viewing angle (except the 9 = 90° case 
which is discussed below). 

An interesting feature is observed at “intermediate” frequencies (~ 6 GHz, 
the second row of Figure [5]), i.e., in the optically thick range, but not far from 
the spectral peak; the optical depth r here exceeds unity, but is not very large 
(t « 10). We can see that the intensity oscillations at these frequencies are 
shifted by about 1/4 wavelength with respect to the oscillations at higher and 
lower frequencies. This is caused by the wave-induced variations of the viewing 
angle: at the intermediate frequencies, these variations (albeit very small) be¬ 
come the main factor modulating the intensity. In contrast, at higher or lower 
frequencies (with r < 1 or t 1), the effect of changing the magnetic field and 
electron concentration dominates. At the same time, the modulation depth at 
the intermediate frequencies reaches its minimum (in comparison with other 
frequency ranges). 

The oscillating magnetic field direction (and hence the oscillating viewing 
angle) also produces some peculiarities in the parts of the loop where the viewing 
angle is close to 90°. Changing the sign of the magnetic field projection on 
the line of sight results in changing the dominating electromagnetic wave mode 
which affects the intensity. As a result, in the optically thick frequency range 
the MHD wave induces intensity variations of opposite signs at different sides 
of the 9 = 90° stripe; the modulation depth increases sharply. In the optically 
thin frequency range the intensity variations at different sides of the 9 = 90° 
stripe have the same sign and the modulation depth demonstrates only minor 
variations. It should be noted, however, that the mentioned peculiarities occur 
only in very narrow stripes and thus their detection in real observations requires 
very high spatial resolution. 

Similar peculiarities (although without significant changes of the modulation 
depth) can be seen for the “Limb A” orientation, in the optically thick frequency 
range, at the distance of about s ~ 2” from the upper footpoint. They occur due 
to the specific loop topology: the optically thick emission is produced near the 
loop boundary and, for the “Limb A” orientation (see Figured]), the dependence 
of the viewing angle at the loop boundary surface on the coordinate s can become 
discontinuous; for other loop orientations, this effect is not observed. 

Figure ini demonstrates the simulation results for the propagating wave. Since 
the structure of the standing and propagating sausage waves is the same, all 
the above conclusions about the radio intensity oscillations remain valid; the 
maximum modulation depths shown in Figuredjare slightly higher than in Figure 
[5| Using Figure [6] we can track propagation of the MHD waves and estimate 
their speed. For the low-density model considered here (see Table [T|), the phase 
speed of the sausage waves equals 8070 km s“^; the visible propagation speed 
of the radio brightenings/darkenings for the “Limb B” orientation is very close 
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to this value. However, due to the varying viewing angle, the projected (visible) 
wave speed varies along the loop; this effect is especially important for the “Disk 
center” and “Limb A” orientations. The amplitude of the propagating radio 
intensity variations also varies along the loop with the most notable effects oc¬ 
curring near the regions of transverse (with respect to the line of sight) magnetic 
field: in these regions, the bright/dark stripes in Figured become discontinuous, 
although the MHD wave itself has no peculiarities; the bright/dark stripes in 
the time-distance plots (in the optically thick frequency range) can also become 
discontinuous in the regions where the viewing angle at the loop boundary 
sharply depends on the coordinates in the image plane. Thus when interpreting 
the radio observations with quasi-periodic pulsations, a special attention should 
be paid to the 0 = 90° layers and/or the regions where the visible loop boundary 
demonstrates topological peculiarities, since they can introduce additional phase 
shifts which do not actually reflect the properties of the underlying MHD waves. 

Assuming that the relative intensity variations are approximately propor¬ 
tional to the relative MHD wave amplitude, we can estimate that the MHD waves 
with the amplitude SB/ (B) of about 40-50% (with other source parameters be¬ 
ing the same as in Figures [Sjni) should produce the intensity modulation depth of 
up to 100%. However, simulating the MHD waves with such amplitudes requires 
taking into account the nonlinear effects. As will be shown in Sections id.4. 211314.31 
the emission modulation depth becomes higher for softer electron beams or for 
higher plasma densities at low frequencies (when the Razin effect is important), 
although this increase of the modulation depth is always accompanied by a sharp 
decrease of the intensity. 

3.4. Parametric Study 

We now consider how the fluctuations of the radio emission depend on various 
parameters of the flaring loop. 

3.4-1. Effect of the Wave Number 

Figure [7] demonstrates the time-distance plots of the intensity variation for the 
models with increased numbers of MHD waves in the magnetic loop (four and 
six waves, which corresponds to the loop lengths of 11300 and 16950 km, re¬ 
spectively); the simulation results shown correspond to the low-density model, 
nonthermal electrons spectral index of <5 = 4, and the optically thin emission (at 
the frequency of 25 GHz). As already shown above, the standing and propagating 
MHD waves affect the radio emission in a qualitatively similar way; therefore 
only the results for the propagating wave are shown. 

Our simulations have revealed no significant effect of the number of the 
MHD waves in the loop, i.e., both the emission modulation depth and the 
phase relations between the emissions at different frequencies remain the same. 
The only difference can be noticed in the regions where the loop axis becomes 
nearly parallel to the line of sight, i.e., near the loop footpoints for the “Disk 
center” orientation and near the loop top for the “Limb A” orientation: since 
the length of the formation region of the optically thin emission can exceed 
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Figure 7. Time-distance plots of the intensity variation at the loop axis (same as in Figures 
[5l6ll at the frequency of 25 GHz for different loop orientations and different numbers of 
wavelengths (four and six) in the loop. The low-density model (propagating wave) is used 
and the electron spectral index is 5 = 4. 


an MHD wavelength, the effects of variation of the magnetic field and electron 
concentration in different parts of the emission source (with opposite phases) can 
compensate each other. As a result, the amplitude of the intensity oscillations 
decreases with a decrease of the relative (with respect to the loop length) MHD 
wavelength. This effect is not observed in the optically thick frequency range. 

It can be seen in Figure [7] (as well as in Figures |4][6]) that the MHD ob¬ 
servations are best resolved when the magnetic tube is observed in a nearly 
perpendicular direction. Indeed, for the “Limb B” orientation we can see one- 
to-one correspondence between the MHD waves and the spatial and temporal 
structure of the radio emission. For other orientations, this correspondence is 
less straightforward and some MHD wave nodes can even be missing in the 
radio emission; as expected, the radio emission cannot accurately reflect the 
MHD wave structure in the regions where the viewing angle with respect to the 
magnetic field 9 is small. While it is not easy to estimate the number of observed 
MHD wave nodes in a general case, we can expect (and this is confirmed by the 
simulations) that the MHD wave nodes become missing in the radio emission 
in the regions with tan0 < R/^, where R is the loop radius and A is the MHD 
wavelength; this happens either due to the averaging along the line of sight (in 
the optically thin frequency range) or due to the strong reabsorption of the 
emission during propagation (in the optically thick frequency range). Note also 
that, as has been discussed above, the 6 = 90° layers can introduce the patterns 
in the radio emission (at low frequencies) similar to additional MHD waves. 
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Figure 8. Average modulation depth of the intensity at the loop axis for different frequencies 
and different spectral indices of the energetic electrons. The low-density model (propagating 
wave) and the “Limb B” orientation are considered. 


3.4-2. Effect of the Spectral Index of Energetic Electrons 


The energy spectrum of the radio-emitting electrons affects greatly the gyrosyn- 
chrotron emission. As can be seen from Figure [3j softer electron beams (with 
larger spectral indices 5) produce less intense emission with steeper spectrum 
slope in the optically thin frequency range; the spectral peak shifts to lower 
frequencies. We have found that variation of the spectral index 8 does not change 
the simulated radio images (like in Figure |4]) qualitatively. However, there are 
quantitative changes in the emission time profiles which may affect detectability 
of the MHD oscillations. 

Figure [S] demonstrates the average modulation depths of the intensity (for 
the propagating MHD wave, three sausage waves per loop, low-density model, 
and “Limb B” orientation); the emission at the loop axis is considered and the 
values are plotted as functions of t he coordinate along the loop. Th e average 
modulation depth e is calculated as ( Mossessian and Fleishmiml . 2012 1 


it is related to the maximum modulation depth A (shown in Figures [5][7]) as 
A ~ e^/2. We can see that in the optically thin frequency range (Figure [H]d) 
an increase in the spectral index results firstly in an increase in the modula¬ 
tion depth. As a result, the oscillation-modulated radio images become more 
contrasted and the periodic variations of the emission become easier to detect. 
For even higher electron spectral indices {8 > 10 for the considered frequency of 
25 GHz), the modulation depth decreases sharply, because the gyrosynchrotron 
emission becomes too weak and the emission is produced mainly by the free-free 
mechanism that is weakly affected by the MHD oscillations. In the optically 
thick frequency range (Figure |B^) the modulation depths are low (< 0.5%) and 
decrease with increasing 8. 

As mentioned above, Figure[H]corresponds to the propagating wave and “Limb 
B” orientation; for other orientations and/or for the standing wave, the above 
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conclusions remain the same. Assuming that the relative intensity variations 
are approximately proportional to the relative MHD wave amplitude, we can 
estimate that for the electron spectral index of i5 = 8 and the MHD wave ampli¬ 
tude bBj {B) of about 25-30%, the average (e) and maximum (A) modulation 
depths in the optically thin frequency range should reach Ri 70% and Ri 100%, 
respectively; the MHD waves with such amplitudes seem to be feasible. 

3 . 4 . 3 . Effect of the Plasma Density 


Finally, we consider the case when the s pectral peak of the gyrosynchrotron radio 
emission is formed by the Razin effect ( Razinl . 1960bl 3) . This effect takes place 
in a relatively high-density plasma and results in a considerable suppression of 
the gyrosynchrotron emission at the frequencies below 


2fl 

3/b’ 


( 6 ) 


where /p and /b are the plasma and cyclotron frequencies in the emission source, 
respectively. At the frequencies below /r, the variations of the thermal plasma 
density should become the main factor affecting the radio intensity. 

As an illustrative example, we use the high-density model (see Table [T|). As 
can be seen from Figure [21 in this case the Razin frequency /r, is about 4 GHz. 
Thus at the frequencies well above 4 GHz the intensity variations should be 
caused mainly by the oscillations of the magnetic field and energetic electron 
concentration (like in the low-density model), at the frequencies below 4 GHz 
the intensity variations should be caused mainly by the Razin effect {i.e., by the 
oscillating plasma density), while at the frequencies around the spectral peak 
(r! 4 GHz) both factors should be important. 

Figure|9]demonstrates the simulation results (time-distance plots of the inten¬ 
sity variation) at three different frequencies for the high-density model, standing 
sausage wave, and the nonthermal electrons spectral index of (5 = 4. In the 
optically thin frequency range, the modulation depth of the intensity is larger 
than that for the low-density model due to the model parameters used (the 
relative MHD wave amplitude for the high-density model is larger, see Tabled]). 
The modulation depth reaches its minimum at a frequency below the spectral 
peak (r 2 GHz, depending on the viewing angle) and then rapidly increases with 
decreasing frequency; however, as can be seen in Figure [S] the intensity at those 
frequencies also becomes very low, so that detection of the intensity oscillations 
in observations would be difficult. 

It can be seen in Figure|9]that for the “Limb B” orientation (when the viewing 
angle is close to 90°) the radio oscillations in the optically thin frequency range 
(> 10 GHz) and near the spectral peak (r 4 GHz) are synchronous (in phase), 
while the low-frequency emission (< 2 GHz) oscillates in anti-phase with the 
emissions at higher frequencies; the phase reversal occurs (for 9 r 90°) at R 2 
GHz. Clearly, this result indicates that at high frequencies (above and even 
slightly below the spectral peak) the intensity variations are caused by the 
variations of the magnetic field and nonthermal electrons concentration, while 
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Figure 9. Time-distance plots of the intensity variation at the loop axis (same as in Figures 
EHZli for the high-density model. Standing wave is considered and the electron spectral index 
is 5 = 4. 


at lower frequencies (< 2 GHz) the Razin effect dominates; these factors affect 
the emission in opposite ways. 

A similar picture can be seen for two other loop orientations in the regions 
where the viewing angle is close to 90° {i.e., near the loop top for the “Disk 
center” orientation and near the footpoints for the “Limb A” orientation). For 
smaller viewing angles the picture is more complicated, since both the modula¬ 
tion depth and the above-mentioned frequency of the oscillation phase reversal 
are angle-dependent: the frequency of the phase reversal increases with a decrease 
of the viewing angle. As a result, if the viewing angle varies gradually along the 
loop, we can see either correlation or anti-correlation between the oscillations at 
two given frequencies, depending on the coordinate _ 

Th e above results contradict the conclusions of Mossessian and Fleishm^ 
( 120121 ) who have found that the oscillations of the high- and low-frequency 
emissions (in the case when the Razin effect is important) should occur strictly 
in phase. As mentioned before (see Section IXiTI) . this is because we consider a 
slightly different model for the variation of the parameters in the sausage wave, 
as well as the spatially i nhomogeneous emission source. Our results also dis agree 
with the conclusions of Reznikova. Antolin. and Van Doorsselaere ( 2014h who 
have found that the oscillations of the high- and low-frequency emissions occur 
mostly in phase (although t his conclusion was dependent on the viewing angle) . 
Most likely, this is because iReznikova. Antolin. and Van Doorsselaera ( 2014l 'l 
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considered the emission with the frequencies below the spectral peak but above 
the “phase reversal frequency” so that the effect of the Razin suppression was 
too weak. 


4. Conclusions 

We have simulated the modulation of the gyrosynchrotron radio emission by the 
sausage-mode MHD oscillations in a semi-circular magnetic loop, in application 
to the solar flares. Although the model used does not include the variation of the 
magnetic field and plasma parameters with height, it has allowed us to study 
the effect of varying the viewing angle along the loop. We have considered three 
loop orientations and different parameters of the magnetic field, plasma, and 
energetic particles. In particular, we have studied in detail the phase relations 
between the spatially resolved intensity oscillations at different frequencies. The 
results can be summarized as: 

• In the model with a relatively low plasma density (when the Razin effect 
is negligible), the high-frequency (optically thin) and low-frequency (op¬ 
tically thick) emissions oscillate, as a rule, in phase. An exception is the 
optically thick emission at the frequencies just below the spectral peak, 
whose oscillations are shifted by « 1/4 MHD wavelength with respect to 
the oscillations at other frequencies. 

• In the model with a relatively high plasma density (when the Razin effect 
is important), the emissions at high and low frequencies oscillate in anti¬ 
phase. The reversal of the oscillation phase occurs, as a rule, at a frequency 
well below the spectral peak; this frequency depends on the viewing angle. 

• In the optically thin frequency range, both the unperturbed intensity and 
the absolute and relative intensity variations are the highest in the regions 
where the local magnetic field is nearly perpendicular to the line of sight. 
In contrast, in the optically thick frequency range all the mentioned values 
are the highest in the regions with small viewing angles. 

• In the optically thin frequency range, the intensity oscillations are more 
pronounced for softer electron beams (ie., with higher spectral indices). In 
the optically thick frequency range the trend is opposite. 

• In the optically thick frequency range, the intensity oscillations can have 
sharp phase shifts between the adjacent regions due to the presence of 
layers of transverse (with respect to the line of sight) magnetic field and 
other topological peculiarities. 

• Averaging over the visible source area and/or along the line of sight re¬ 
duces the amplitude of the intensity oscillations; this effect becomes more 
important with an increase in the number of the MHD harmonics in the 
loop. 

The latter result means that, due to a limited spatial resolution, the real observa¬ 
tions can differ from the ideal simulations presented here. Nevertheless, we expect 
that if an instrument spatial resolution is better than a MHD wavelength, this 
will allow us to detect (or rule out) the above-mentioned phase relations between 
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the oscillations at different frequencies; another requirement is the ability to 
perform imaging observations at many frequencies. We anticipate that the forth¬ 
coming multiwavelength radioheliographs (such as the Expanded Owens Valley 
Solar Array, Upgraded Siberian Solar Radio Telescope, and Chine s e Spe ctral 
Radioheliograph; see the review of Nakariakov. Kashanova. and YanI ( 2014l l and 
references therein) will be able to provide the necessary data. 
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Appendix 

A. Toroidal Coordinates 


As was mentioned above, we consider the MHD waves in an overdense cylinder; 
the system can be described by the cylindrical coordinates: the distance along 
the cylinder axis 2 , the distance from this axis r, and the azimuthal angle ip (see 
Figure [TUhl. We assume that the cylinder has the length L {i.e., 0 < z < L) 
and the radius R (the radio emission is produced only inside the cylinder, i.e., 
at r < i?). 

The cylinder is assumed to be bent into a semi-circular (semi-toroidal) loop 
(see Figure fTHb l: the loop curvature radius is H = L/tt (evidently, we need 
H > R). We assume that initially the loop is oriented as shown in Figure fTOb. 
i.e., it is vertical, symmetric with respect to the x'z' and y'z' planes, and its 
axial line lies in the x'z' plane. Then the cylindrical coordinates (z, r, tp) are 
related to the Cartesian coordinates {x',y',z') as 


and 


x' = (rcosv?-I-i?) cosy, (7) 

y' = rsintp, (8) 

z' = {rcosp + H)smx, (9) 


tany 


tan p 


{x' - H cosxf + {y'f + {z' - H sinx)2, 

_ ][_ _ 

x' cos X -I- z' sin X — ’ 


( 10 ) 

( 11 ) 

( 12 ) 


where x = z/L. 

Assume that a vector A at a given point is defined by its three components 
[Az, Ar, A,p) in the cylindrical coordinate system. After transformation to the 
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Figure 10. a) Sketch of the model plasma cylinder, b) Sketch of the model semi-toroidal loop. 


semi-circular loop, the vector can be expressed as A' = -|- A,^a^, 

where the unit vectors r the Cartesian coordinates are given by 


(-sinx,0,cosx), 

(13) 

(cos ip cos X, sin tp, cos (/? sin x), 

(14) 

(— sin ip cos X, cos ip, — sin p sin x). 

(15) 


The above formulae are valid in a general case. In this work we consider the az- 
imuthally symmetric case when all parameters are independent of the azimuthal 
angle and the azimuthal components of all vectors are equal to zero. 

The orientation shown in Figure [TUb corresponds to the vertical loop located 
at the center of the solar disk, with the axis z' being the line of sight and the 
plane x'z' being the equatorial plane. All other orientations (see, e.g., Figure [2]) 
are obtained by rotating the loop firstly by the angle ip around the local vertical 
(z' axis) and then by the angle A around the solar rotation axis. 
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